A simple model of Coulomb disorder and screening in graphene 
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We suggest a simple model of disorder in graphene assuming that there are randomly distributed 
positive and negative centers with equal concentration N/2 in the bulk of silicon oxide substrate. We 
show that at zero gate voltage such disorder creates the two-dimensional concentrations no ~ N 2 ^ 
of electrons and holes in the graphene sample. Electrons and holes reside in alternating in space 
puddles of the size _Ro ~ iV -1 / 3 . A typical puddle has only one or two carriers in qualitative 
agreement with the recent scanning single electron transistor experiment. 



Recent experiments on the gate voltage dependent 2D 



1.2.3 



at- 



transport of graphene deposited at SiO substrate^ 
tracted almost unprecedented theoretical attention (see, 
for example2i£i£iL£i2, and references therein). Clean 
graphene has zero mass Dirac relativistic spectrum near 
the Fermi level (e = 0), where energies e of both electrons 
and holes behave as 

£ = ±hvk. (1) 

Here the velocity v ~ 10 6 m/s and k is absolute value of 
electron wave vector. Density of states for such a (four 
times degenerate) spectrum v(e) — 2e/-Kh 2 v 2 vanishes at 
e = 0. Only states with a negative energy are filled in a 
clean graphene sample. 

Experimental measurements of the electrical conduc- 
tivity and the Hall effect voltage strongly indicate pres- 
ence of disorder. Most of the theoretical discussion is fo- 
cused around observed phenomenon of saturation of the 
linear dependence of two-dimensional conductivity a on 
the gate voltage at small gate voltages, where the con- 
ductivity reaches the minimum value a = a m i n = Ce 2 /h. 
Particular attention is attracted by the constant C in this 
formula^. 

In this note we concentrate on the width of the con- 
ductivity minimum. Using the linear relation between 
gate voltage and concentration of electrons we charac- 
terize this width by the concentration n = hq, where 
saturation starts. The hint to the meaning of n = hq fol- 
lows from the Hall effect data, well fitted by the so-called 
two band model, which assumes simultaneous presence 
of electrons and holes even at zero gate voltage^. It was 
suggested^ii^ that one can understand these phenom- 
ena assuming that potential of charged impurities moves 
the Dirac point up and down in different points of space 
creating alternating in space electron and hole puddles. 
Theoretical self-consistent calculations of the random po- 
tential and carrier distributional^ so far were based on 
models of two-dimensional distribution of charged impu- 
rities, which require at least two parameters, the total 
two-dimensional concentration of impurities and width 
of the layer near the surface they reside. 

Here we would like to suggest a simple model of distri- 
bution of charged impurities and a theory of their non- 
linear screening. Below we assume that there is three- 
dimensional concentration N/2 of both positive and neg- 
ative centers randomly distributed in the bulk of SiO 



substrate and estimate all mentioned above quantities in 
terms of the single parameter of the model N. We esti- 
mate the concentrations of electrons and holes no at the 
zero gate voltage (or, more exactly, in the center of mini- 
mum of conductivity) , the characteristic spatial scale Rq 
of nonlinearly screened potential (which gives the size 
of the typical electron or hole puddle and the distance 
between neighboring electron puddles) and the typical 
number of carriers in a puddle in terms of N. For this 
purpose we study screening of the random potential of 
3D impurities by a graphene sample. 

Nonlinear screening of the random potential of charges 
impurities by electron and hole puddles was first studied 
for three-dimensional totally compensated semiconduc- 
tor s 11 ! 12 ' 13 . Following these works we start from calcu- 
lation of the random potential with the spatial scale R. 
We imagine that all oxide is divided in cubes with the 
edge length R. Each cube has average number of im- 
purities NR 3 and fluctuating excessive charge of either 
sign e(NR 3 ) 1 / 2 . This charge creates the random poten- 
tial fluctuating from a cube to cube with the amplitude 
of the order of 

V(R) = {e/nR^NR 3 ) 1 ' 2 = {e/K,)(NR)^ 2 , (2) 

where K ~ 2.5 is the effective dielectric constant at the 
oxide surface. Here and everywhere below we drop nu- 
merical coefficients making only order of magnitude es- 
timates. All the cubes adjacent to the surface of oxide 
apply this potential to the graphene sample. At zero gate 
voltage (in conductivity minimum) there are no electrons 
and holes without this potential. The random potential 
moves e = point up and down and creates electrons 
in its wells and holes on its hills. Let us concentrate on 
electrons. Electrons screen a typical well of the size R al- 
most completely if they can fit in this well in the number 
equal to the well charge. The concentration of electrons 
necessary for neutralization of the well charge is 



n{R) = {NR 3 yl 2 /R 2 = {N/R) 1 ' 2 . 



(3) 



Let us assume that the wavelength of these electrons 
n" 1 / 2 is smaller than R, or nR 2 ;§> 1, so that the Fermi 
energy of them electrons can be estimated in the Thomas- 
Fermi approximation, 



E F (R) = Hvn 1 / 2 = (N/R) 1 / 4 . 



(4) 
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If Ep(R) <C eV(R) all electrons compensating the well 
charge are localized inside the well close to its bottom 
and are fragmented in space by the random potential 
of even smaller sizes. Thus we are dealing with nonlin- 
ear screening when concentration of electrons is strongly 
nonuniform. 

The process of fragmentation of electron liquid contin- 
ues until eV(R) 3> Ep(R). It stops at such R = Rq, 
where Ef(Rq) = eV(Ro). This gives equation for R 

hv(N/R ) 1/4 = (e 2 / K )(NR ) 1/2 , (5) 

which yields Rq = N-^ 3 (e 2 /kHv)- 4 / 3 . Using v = 10 8 
cm/s we get that e 2 /nhv ~ 0.8 and R ~ A" 1 / 3 . The 
smallest wells of the size Rq are not fragmented and form 
what one may call electron puddles. A typical puddle has 
one or two electrons. Groups of one of two negative impu- 
rities at the distance iV -1 / 3 from graphene form hole pud- 
dles. Condition nR 2 S> 1 used above is only marginally 
correct, so that the whole calculation is only an order of 
magnitude estimate. 

Now we can estimate the two-dimensional concentra- 
tions of electrons and holes. Each of them occupy half of 
space with the concentration 

n ~ n(R ) ~ N 2 ' 3 . (6) 

Eq. ([fj]) is the main result of our paper. It gives the two- 
dimensional concentration of coexisting electrons and 
holes no in graphene in terms of a single parameter of 
our model, three-dimensional concentration N of charged 
impurities. It could be easily anticipated because this is 
the only formula with necessary dimensionality one can 
make from N . There are no large or small dimensionless 
parameters in the system, which could play the role of 
large or small coefficient in Eq. ©. Indeed, the single 
dimensionless parameter e 2 / ' kTiv is close to unity. 

Above we are using the term puddle size for the small- 
est scale of fragmentation of electron and hole density. 
Although electrons and holes are really located in such 
small puddles there are fluctuations of the electron den- 
sity at larger scales R 3> N~ x / 3 . At such scales the elec- 
tron density follows Gaussian fluctuations of the density 
of impurity charges. For example, for domains of size R 
the electron concentration n(R) is given by Eq. © and 
is much smaller than no- Although some people may re- 
fer to these domains also as "puddles", they are actually 
small gaussian fluctuations of the density of puddles of 
one or two electrons we are talking about. In other words, 
in the plane of graphene large scales of the random po- 
tential R 3> Rq are screened linearly, by small variations 
n(R) around the average concentration no of electrons 
and holes situated in puddles. One can check that at 
e 2 /kHv ~ 1 the linear screening radius corresponding to 
the two-dimensional concentration no is equal to Rq. 

Recently both the concentration of carriers no and the 
characteristic size of puddles Rq , were estimated by com- 



pressibility measurements based on the use of the scan- 
ning single electron transistor—. The authors arrived at 
no = 2.3 10 11 cm -2 and Rq = 30 nm. These numbers 
result approximately in uoRq ~ 2 electrons per puddle 
in a good agreement with our estimates. Using Eq. ([6]) 
we can also estimate that the three-dimensional concen- 
tration of charges in the silicon oxide substrate N ~ 10 17 
cm -3 , what is not unreasonable. 

We would like to emphasize that while working with 
limited resolution authors of Ref, — saw domains of fluc- 
tuations of density of electrons with larger sizes R Rq, 
but with much smaller concentrations n(R) <C no- They 
noticed that n(R) increases with the decreasing scale R, 
and improving experimental resolution traced observed 
gaussian dependence of n(R) to the smallest scale cited 
above. Their observations are in agreement with our pic- 
ture. 

In the recent paper— authors suggested a theory of 
graphene conductivity assuming existence of macroscopic 
puddles totally filled by electrons and holes and separated 
by narrow p-n junctions. In other words, they assumed 
that each electron puddle has many electrons. We see 
that in our model macroscopic approach to puddles is 
not justified. This may lead to some numerical changes 
in calculation of constant C in the expression for a m in 
above and in the shape of the conductivity versus gate 
voltage curve near the minimum, but definitely will not 
change the qualitative conclusions of Refi*. 

Note that fragmentation of the electron density in 
graphene into puddles of one or two electron is not an 
exclusive feature of the model of three-dimensional dis- 
tribution of impurities. Similar results follow from ran- 
dom distribution of charges in a two-dimensional plane 
close to graphene^i 7 -. Note also a difficulty with appli- 
cation of the model of this paper for evaluation of the 
Coulomb scattering limited mobility of electrons at large 
concentrations Ti>no. Indeed at large Fermi wave vec- 
tors ~ n 1 / 2 electrons are strongly scattered only by the 
the Coulomb centers at the distance of the electron wave- 
length n^ 1 / 2 <C N~ x / 3 from the graphene plane. In other 
words, electrons do not "see" more distant centers. As 
a result the mobility should increase proportionally to 
n 1 / 2 . This prediction is in contradiction with the ob- 
served independent on n mobility, which in turn is well 
explained by Coulomb impurities residing in graphene 8 . 
This may be a strong argument for the dominating role 
of charged centers located in graphene in the currently 
studied samples. But even if this is proven, in future 
charges in graphene may be neutralized or eliminated. 
Then one will arrive to the case the bulk charged impu- 
rities discussed here. 

I am grateful to V. Galitski, D. Novikov and K. 
Novoselov for very useful discussions. 
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